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ABSTRACT 


Hcrizontal and vertical thermal structures were examined 
in a region of the Northeastern Pacific Ocean off Northern 
California. The observations were acguired on three cruises 
as fart of the ONK-sponsored Ocean Prediction Through 
Cbservation, Modeling and Analysis (OPTOMA) research 
program, centered ina region of the California Current 
System, ca. 37 to 399N, 124 to 126°W, during June and July, 
1383. 

The horizontal temperature correlation scale was between 
30 and 50km, which was a significant factor when comparisons 
were made between measured horizontal thermal structures and 
those retrieved from Fleet Numerical Oceanography Center's 
{(FNOC) analyses, which had grid spacing of 320km ({TEOTS) and 
40km (EOTS). 

Operational (modeled) analysis fields were not in agree- 
ment with the observed fields. The major discrepancies 
cccurred in the magnitude or the mixed layer depth and the 
Shape of the horizontal temperature fields (maps). 

The operational significance of differences between 
modeled and observed thermal structures was assessed in 
terms of their effect on low-frequency; i.e.,less than 1Khz, 
acoustic propagaticn Ut2 leZang the FACT9H and PE 
Transmission Loss nodéels. 

Mixed layer depth differences produced significant disa- 
greements between direct ranges computed from model and 
observed temperature profiles. The effect was most 
pronounced at higher frequencies and when both source and 
receiver were shallow; i.e., both at 20m. 

A comparison was made between average depth/temperature 
profiles from July, 1982, July, 1983, and FNOC climatology 


to obtain a measure of the effect of interannual variability 
in the domain. This comparison showed that a significant 
temperature anomaly existed inthe upper 400m in 1983 
compared to 1982 due to El Nino and that this anomaly was 
not represented by tke FNOC climatology. 

The differences tetween modeled and measured thermal 
structures are believed to be related to thermal structure 
model resolution, model sensitivity to input data, short 
scales of spatial variability and non-representative clima- 


tology for the domain. 
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I. INTRODUCTION 


A. BACKGROUND AND OBJECTIVES 


The dependence of acoustic propagation on the ccean 
thermal structure and the increased demand from the Navy for 
accurate acoustic forecasts to facilitate planning and 
economical use of resources has led to an increasing 
requirement for a reliable and accurate model to represent 
ocean thermal structure. Such a reguirement iS not easily 
met, however, because of the limited number of observaticns 
in the typical ocean thermal structure data base, approxi- 
mately 150 XBT measurements per day in the northern heni- 
sphere [Clancy et al.,1983]. In view of the upper ocean's 
Spatial and temporal scales of variability such a small data 
bkase represents a Significant handicap to the performance of 
any ocean analysis ncdel. 

The availability of an ocean thermal structure data base 
from a relatively densely sampled area of the Northeastern 
Pacific Ocean (approximately 200km west of point Arena, 
Northern California) made it possible to evaluate the accu- 
racy of the output of two operational ocean thermal struc- 
ture analysis systems. 

The cbjectives of this research are firstly to obtain a 
description of the characteristics of the ocean therzal 
structure from a pcrtion of the ocean (ca.379N to 39°N; 
1249W to 126°W) sampled during the three cruises comprising 
CPTOMAS which tock flace in June and July 1983. (OPTOMA - 
Ocean Prediction Through Observations, Modeling, and 
Analysis, is a joint Harvard/NPS project tintended to 
acquire field data to characterise synoptic scale eddies 


over a domain in the Caiifornia Current off Northern 
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California, and to ‘set-up an eddy-resolving, statistical/ 
dynamical, limited domain, open boundary numerical ocean 
predicticn model.'Tke '5* indicates the fifth series of 
cruises of the project). Secondly, the analysed ocean 
thermal structure is to be used as a ‘sea truth" by which to 
evaluate the accuracy of two real-time ocean thermal anal- 
ysis systems presently in use at the Fleet Numerical 
Cceanography Center (FNOC), namely, the Expanded Ocean 
Thermal Structure (ECTS) analysis andthe Thermodynamic 
Ocean Prediction System (TOPS)-Coupled EOTS (TEOTS) anal- 
ysis. The final objective is to determine the operational 
Significance of deviations in the modeled thermal structure 
from the ‘sea truth by utilizing the output from opera- 
tional acoustic models to provide a guantitative measure of 
the effect of these thermal structure differences on 
acoustic propagation in the ocean. Two low-frequency, 
passive-acoustic models were used to provide a measure of 
the variability of acoustic energy propagation which could 
be expected from the acoustic models alone. These then, 
would provide a benchmark by which the variability due to 
differences between modeled and analysed thermal structures 
could ke judged significant or not. The two acoustic models 
used were the FACT9H (Fast Asymptotic Ccherrent 
Transmission) transmission loss model and the AESD (Acoustic 
Enviromental Support Detachment) PE (Parabolic Eguaticn) 
model [{ Brock,1978 }. 


Be. STUDY AREA AND ITS OCEANOGRAPHY 


The area investigated is located in the Califcrnia 
Current System (CCS), south of the Mendocino Escarpment. The 
regicn 1S approximately bounded by 379 to 39° N and 124° to 
126°W. The bathymetry, Fig. 1.1, Shows the area to be 
offshore of the continental slope, with an average water 


depth of 4000nm. 
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The CCS is composed of four currents: the Califernia 
Current, the Davidson Current, the California Undercurrent, 
and the Southern California Current. Of these four currents 
the southward flowing California Current is the major influ- 
ence in the study dcmain. This current forms the eastern 
boundary current of the large anticyclonic gyre centred near 
the Hawiian Islands. As an extension of the Westwind Drift, 
the California Current occurs between the North Pacific 
atmospheric high pressure system and the semi-permanent 
thermal low positioned over Central California. Contrary to 
the classical picture of a broad, slow, shallow climatolo- 
gical mean eastern boundary current, the CCS is known to 
comprise mesoscale eddies, meandering jets and turbulent 
filaments [Mooers and Robinson, 1984] together with their 
associated complex hcerizontal and vertical thermal struc- 
tures. The mean current is broad (approx. 1000km wide), 
shallow (less than 500m deep) and sluggish moving (approx 
O.9kmyhr). The southward to southeastward winds, which 
prevail during the spring to the fall months south of 40° N, 
produce upwelling aicng the coast, extending 50 to 100km 
offshore. Upwelling begins in February or March and 
continues through August in the study domain. The cold, 
usually high salinity water, brought into the surface layers 
by this frocess, is then moved offshore by the smaller scale 
perterbations in theccs, further complicating the ocean 
thermal structure. The most intense upwelling sites cccur to 
the south of capes and headlands [Reid et al.,1958] with the 
most conspicuous centres of upwelling adjacent to the study 
domain at 359 N and 41° N [Sverdrup et al., 1942 and Reid et 
ley Sele 

Mean monthly upwelling velocities for the west coast of 
North America(Fig 1.2; from Heath (1983)), based on Nelson's 
(1977) wind stress curl data, clearly show upwelling to be 


at its most intense in July off Northern California. Hence, 
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the maximum horizontal and, because of the differing water 
masses involved, vertical variations in thermal structure 
Should be expected to occur in the study domain during the 
period under investigation. 

In addition to the seasonal variations in thermal struc- 
ture, inter-annual variability is also evident. Significant 
increases in temperature off Central and Northern California 
occurred during the winters of 1972, 1976, 1979 and 1983, 
corresponding to El Nino episodes [Breaker, 1983). Of 
particular interest are the effects of the major El Nino 
warming of 1982-83. A brief analysis of the possible opera- 
tional ccnseguences of such an anomolous event is presented 


in Chapter V. 
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Ti. DESCRIPTION OF COMPUTER NODEL 


A. OCEAN ANALYSIS MCDELS 


t. Expanded Ocean Thermal Structure (£EOTS) Analysis 


The Expanded Ccean Thermal Structure (EOTS) analysis 
[Mendenhall et al., 1978] has been the US Navy's operational 
thermal analysis system for the past several years and 
provides ocean thermal structure input to most of the 
acoustic predictions generated at FNOC. 

EOTS is run cn a 63x63 grid for both hemispheric and 
regional coverage. The grid spacing for the hemispheric 
product is 200 to 400km, while that of the regional version 
is typically 40 to 80km. For the SOCAL region, in which the 
study domain resides, the spacings are 320km and 40kn 
respectively. 

The EOTS analysis systen is based on a 
Fields-by-Information-Blending (FIB) methodology [Holl and 
Mendenhall, 1971; Hcll, Cuming and Mendenhall, 1979) and is 
a numerical blending technigue whereby the difference 
rketween a first guess field (in this case the previous anal- 
ysis weighted towards climatology at a specied rate) and the 
assembly field, are eliminated in a weighted, least mean 
Squares sense. The first guess field provides shape infornma- 
tion for the horizontal blending procedure. The assembly 
field, into which all new observations are assembled, is 
obtained in a two-step process. First, the assembly field 
from the previous analysis is trended towards climatology in 
the same manner as tke first guess field, and a weight field 
associated with it is decayed an arbitrary amount. Then, a 
weighted average of the assembly and first guess fields is 
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formed (using a low, specified weight for the first guess 
field) to provide a new assembly field. New observations are 
interpolated to the nearest analysis grid point using shape 
information from the first guess field. At each grid point, 
a weighted average of the interpolated observations and the 
assembly field is formed and a final assembly field deter- 
mined. It iS upon this final assembly field and the first 
guess field that the FIB analysis is performed. The proce-~ 
dure is designed to spfread information from data-rich points 
to data-sparse points and produces a smoothed, horizontally 
and vertically blended field. During the blending, both 
first guess and assemtly fields are weighted towards clima- 
tology, and, in the absence of inputs, they will revert to 
climatclcgy over a period of approximately fifteen days 
[Pollack, personal ccmmunication}. 

The analysis is performed on twenty-six parameters 
at fixed and floating levels from the surface to 400n, Fig 
2.1. Parameter 11s the Primary Layer Depth (PLD) which in 
general equates to the depth of the surface mixed layer. The 
remaining parameters are temperature at selected depths and 
vertical temperature gradients. Parameters 2 to 8 are anal- 
ysed at floating levels defined relative to the PLD. 
Parameters 9 to 26 are associated with fixed levels. Below 
4GOOm the thermal field is constructed from climatological 
archives modified to blend smoothly with the analysed temp- 
erature profiles above. 

Selection of the PLD (point 1) 1s accomplished by 
applying an algorithm to all BT observations made during a 
specified number of days (typically five) prior to the anal- 
ysis time to identify PLD candidates (1i.e., the depth(s) at 
which the criteria prescribed by the PLD algorithm are 
satisfied) from each frofile. The algorithm assigns the PLD 
to the depth(s) (below the surface) of maximum rate of 


curvature in the vertical profile. Then a preliminary 
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three-cycle FIB is performed using the selected FLD 
candidates as input, and the PLD and assembly fields fron 
the previous analysis, trended towards climatology, as the 
first guess and initial assembly fields. Any recent BT's are 
then examined and the PLD candidate from each that is 
Closest to the PLD prceduced by the preliminary FIB analysis 
is selected and entered into a second three-cycle fFIB 


analysis. This produces the final analysed PLD. 


2. TOPS-Coupled EFOTS (TEOTS) Analysis 


The TEOTS system iS an attempt to incorporate 
oceanic fhysics and air-sea interactive processes inte the 
analysis by coupling the ocean thermal analysis to the 
atmospheric forcing via the physics of the TOPS mixed layer 
model. 

The TOPS mixed layer model is the Navyts vertical 
ocean thermal structure forecast system. Given that the 
ocean 1S primarily atmospherically forced, and that mixed 
layer depths have proved *thighly predictable with a variety 
of models't, such as the Denman (1973) or Mellor and Durbin 
(1975) models, TOPS was developed by NORDA (Naval Ccean 
Research and Development Activity); it is based largely on 
the use of rather well-defined atmospheric variables to 
drive a model to describe the relatively data sparse oceans. 
Surface wind, solar heat fluxes, and precipitation fields 
are among those necessary to drive the TOPS/TEOTS models. 
The surface fluxes are supplied by NOGAPS (Naval Operational 
Global Atmospheric Prediction System) and the physics of the 
model are based on the Mellor and Yamada (1974) level-2 
turbulence model with the added effects of advection by 
instantaneous wind drift and climatologically averaged geos- 
trophic currents. The TOPS model produces a twenty-four 
hour forecast for temperature fields at seventeen depths 


from the surface to 500m on a horizontal grid spacing 
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identical to that used in the teots analysis to which it is 
coupled. The main purpose behind the TOPS/TEOTS apfroach is 
to de-couple the objective analysis component (EOTS) SST and 
MLD fields from climatology, and to provide realistic lower 
boundary conditions fcr the atmospheric models. 

TEOTS is coupled to TOPS in cyclical fashion, 
providing initial conditions for a 24hr TOPS forecast that 
is sukseguently fed kack into EOTS asa first guess field 
for the following day's analysis. With the exception of the 
lock-step procedure described above and a different 
prescription for certain parameters, the TEOTS analysis 
system is the same as EOTS, relying on the FIB technique to 
combine the various types of data. Clancy and Pollack (1983) 
give a mcre rigorous description of both EOTS and TEOTS. 

For the study domain EOTS was used for the regional 


(NOCAL) analysis and TEOTS for the hemisphericai analysis. 


B. ACOUSTIC MODELS 
1. The FACTOH Ncdel 


The FACT9H M@Iransmission Loss Model is the Navy 
Interim Standard Transmission Model for ocean regions which 
may be treated with a single sound-speed profile and a flat 
ocean bottom [ Jacobs, 1982). It 1S a fray-acoustics model 
designed for the comfutation of transmission loss as a func- 
tion cf range and freguency for a source and receiver at 
fixed depths. The classical ray treatment has been augmented 
with the addition of higher order asymptotic corrections in 
the vicinity of caustics, and the phase addition of certain 
ray paths. 

Wave theory for acoustic propagation states that 
there 1s a frequency below which no 'trapped' nornal mode 
may exist for a given channel or duct (waveguide) [Cofpens 
et al., 1980]. The FACT9H model attempts to translate this 
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condition to ray-acoustics by applying an algorithm [ Clay, 
1968] to characterize the gross features of ducted propaga- 
ti@m. The intensity in the surface duct is found from the 
principle of conservation of energy modified by additional 
losses (proportional to range) caused by duct leakage and 
rough-surface scattering of energy from the duct [Marsh and 
Schulkin, 1967]. For both source and receiver within the 
surface duct the transmission loss (TL) asa _ function of 


range (R) is given by Eqn. 2.1 
TL(R)=22 + 10log(R.@) + D.R (eqn 2.1) 


where @ is the angle for the ray at the surface which just 
grazes the bottom of the duct and b contains the duct- 


leakage and rough-surface losses,Egn. 2.2 
b=14.88 x 108 £-38 g-¥3 z-3 + W £¥2 z-¥2 (eqn 2.2) 


In the above expressicn f 1s the frequency, g the magnitude 
of the below layer gradient, 2zthe duct depth, andwWisa 
factor determined by surface wave height. 

The 9H is a model designator and reflects the stage 
of development of the model physics (FACT10 is presently in 
operation at FNOC). 


2. The PE Model 


The PE Model replaces the reduced elliptic wave 
equation with a parakclic partial differential equation that 
can be integrated numerically using the Tappert-Hardin 
Split-step Fourier Algorithm [ Brock, 1978). The parabolic 
wave equation includes diffraction and all other full-wave 
effects as well as range depedent environments. The entire 
range and depth dependent acoustic field is computed as the 


soluticn is marched forward in range. The model assumes a 
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flat pressure release ocean surface and a vanishing field at 
the depth of the finite Fourier transform (i.e., for bottom 
grazing angles greater than the maximum permissible 339, the 
acoustic energy is smoothly attenuated to zero and. back 
reflection of acoustic energy into the water column is 
avoided). For grazing angles less than the maxinun, a 
pseudo radiation condition is introduced at the water-bottonm 
interface by smoothly attenuating the field. Since the error 
in the parabolic approximation increases as angles increase 
from the horizontal, steep bottom slopes can cause inaccura- 
cies. The model is primarily considered useful for 
predicting low-freguency acoustic propagation of energy 
along waterborne or shallow-angle, bottom-bounce paths in 
deep water. 
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III. DATA BASE 


A. DATA ACQUISITION 


The XBT and CTD data from the surface to 500m acquired 
from the OPTOMAS cruises formed the primary source of field 
data. CPICMAS consisted of three separate cruises undertaken 
ketween 15 June and 19 July, 1983 by R/ V ACANIA. The three 
cruises, designated AI, AII and AIII each sampled approxi- 
mately the same area of the CCS, Fig. 3.1 

A secondary source of field data was provided by deep 
CTD casts from the cruise undertaken by the R/V DE STEIGUER 
retween 10 and 21 June, 1983. A section of this cruise 
Sampled the OPTOMAS area and provided CTD data to 3000m. 
Satellite imagery was also used, when available, to provide 
a gualitative measure of the sea surface temperature 
distribution. 

Temperature profiles, to be compared later with the 
field data, were retrieved from the archives of FNOC for the 
Same area and period. Both EOTS and TEOTS depth/temperature 
profiles for the beginning and end of each OPTOMAS cruise 
were obtained for selected positions, Fig. 3.2 . 


Details of the data base are given in Table I . 


Be. DESCRIPTION OF DATA AND DATA PROCESSING 


1. ACANIA Data 


ALGO IG APRS 


The positions of the XBT/CTD casts for the Acania 
cruises are documented, Figs. 3.3 to 3.5. All XBT/CTD Gata 
were edited thus eliminating any obviously erroneous 
profiles. A secondary visual inspection of each profile was 
undertaken with the result that a final data base for each 


cruise was obtained, Table I 
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TABLE I 


| 
Details of Data Base 


SOURCE PERIOD NO. OF ACCEPTABLE RANGE OF 
PROFILES PROFILES 
OPTOMA AI 6/ 18-6/22 74 0-500m 
OFICMA AII 6/29-7/4 89 0-500mn 
OPTOMA AIII 7/13-7/19 86 0-500n 
DE STEIGUER 6/10-6/21 8 0-3000m 
6/ 16,6/23, 
EOTS/TEOTS 6/30,7/7 42 0-40 Om 
T7105 772 


The XBT's were then interpolated in the vertical to 
Standard 5m depths from the surface to 500m. The loss of 
temperature microstructure resulting from this interpolation 
was assumed unimportant in view of the goals of this study. 
These e€ditec, ‘smoothed" depth/temperature profiles fcrmed 
the field data upon which the later analysis was to be 
based. 

Information from each cruise waS assumed ‘quasi- 
Synoptic (i.e. the region was believed to have varied some- 
What over the period of sampling). This assumption is 
weakest in the near surface layers where the time scale of 
variability can be of the order of a day or 1éess. The 
validity of this assumption in the context of the present 
study is addressed in Chapter V. 


2- DE STEIGUER Data 


The DE STEIGUFR sampling provided CTD casts from the 


Surface to a maximum depth of 3000m. These profiles were 
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also interpolated vertically to 5m depth intervals and, 
based onthe assumption that the water column was near 
stationary below the thermocline (approximately 500m), they 
were used to extend the ACANIA and FNOC profiles to a depth 
of 3000m. These profiles were then extrapolated linearly to 
the ocean bottom (assumed flat at 4000m). 

An average salinity/depth profile was calculated 
from the ‘in-area‘' [TE STEIGUER CTD casts and this average 
profile was assumed as the standard salinity profile for all 
sound speed calculaticns. 


3. Satellite Imagery 


Advanced Very High-Resolution Radiometer (AVHRR) 
infrared images from the NOAA-7 satellite were used to gqual- 
itatively show the mrcre obvious features of the SST field. 
These images were obtained on an ‘aS available* basis during 
each of the ACANIA cruises. Details of the major features 
Shown in these images were sketched, Figs. 3.6 to 3.8 The 
Satellite images were also used to qualitively assess the 


assumption of guaSi-synopticity. 


4. ENOC Archived Data 


Temperature profiles for the standard positions, 
Fig. 3.2 , were extracted from the archives of FNOC. ‘These 
profiles were used fcr both EOTS and TEOTS analyses on the 
dates closest to the commencement and termination of each 
ACANIA cruise. 
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A. JTCATA PROCESSING 


To facilitate ccm&parison with FNOC's data fields, the 
study demain was divided into a 7x7, Sguare grid witha 
grid-spacing of 20km and the irregularly spaced xXBT/CTD 
profiles interpolated onto this grid. 

The horizontal ccrrelation of the temperature field at 
various depth intervals in the top 150m of the water column 
were evaluated. The correlation at various levels displayed 
Similar chracteristics showing zero crossings occurring at 
approximately 80km with a0O.5 correlation at approximately 
30km (approximately the e-folding distance between the 
origin and the zero crossing distance). It was, therefore, 
reasonable to assume that the temperature variations at 
points within 30 km cf a given position within the domain 
were related (i.e. all points within 30km could be expected 
to influence a given fyoint). As an example, the SST hori- 
zontal ccrrelations for each cruise are presented, Figs. 4.1 
to 4.3 . The correlation information thus derived was used 
when interpolating the field data onto the grid. Only 
observations within 30km of each grid point were used and 
each observation value was weighted according to the inverse 
of its distance from the grid point. The interpolated data 
thus fcrmed were then automatically contoured at FNOC temp- 
erature levels, Fig. 2.1 The mixed layer depth (MLD) calcu- 
lated frem the field data was equated to the PLD in the FNOC 
analyses (point 1 Fig. 2.1) and was defined as the shal- 
lowest depth below which a negative temperature gradient was 
found with a magnitude greater than 39°C /100m over’ two 
consecutive 5m depth intervals. The algorithm used to 
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calculate the MLD was designed to yeild the true MLD even in 
the presence of strong, near-surface negative temperature 
gradients caused by diurnal heating the “afternoon effect" 
which seldom exceeded 7m in depth. MLD's calculated by this 
method agreed well with manually determined MLD's for the 
sampled profiles. FNOC EOTS and TEOTS profiles were also 
interpolated and contoured in a similar manner and on the 
same grid. The MLI computed for the model profiles using 
the algorithm described above was in agreement with the PID. 

By contouring the fields in the manner described, it was 
possible to extract a depth/temperature profile from the 
surface to 4000m for any position within the study domain. 

Throughout the study domain, temperature variaticns at 
GOOm were small enough (<0.5°9C) for all profiles (including 
those frovided by FNOC) to be smoothly blended be€etween 
shallow and deep profiles above the depth of the deep sound 
channel axis (ca. 500m). 

The feasibility of using discrete depth/temperature 
levels, with depth separations dictated by the EOTS analysis 
levels, to represent the vertical temperature structure of 
the area was tested Ly examining the vertical pattern corre- 
lation between temperatures in the near surface layers (the 
zone of Maximum variation of temperature with depth), Figs. 
4.4 to 4.6 . These pattern correlations are similar for each 
cruise and have zerc crossings at depths of the order of 
300m when the SST is correlated with deeper levels. Such a 
pattern correlation gives an good indication of the frob- 
ability that variations in temperature at one depth are 
related to variations in temperature at another depth. At 
the e-folding depth (approximately 100m), the correlation is 
approximately 0.7. It was therefore considered reasonable to 
represent the vertical temperature structure at the discrete 
levels dictated by the EOTS analysis. The EOTS depth 


interval between temperature levels in the near surface 
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Figure 4.1 SST Horizontal Correlation-Cruise AI. 
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Figure 4.3 SST Horizontal Correlation-Cruise AIII. 
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layer is of the order of 50m which yields a pattern correla- 
tion coefficient ranging from 0.6 for cruise AII to 0.8 for 
cruises AI and AIII. Below the mixed layer (approximately 
50m ) the correlaticn improves to 0.8 or better between 
adjacent contoured levels. 

The smoothing effect of deriving profiles as described 
above (i.e. the higher frequency fluctuations of temperature 
with depth are removed) was not considered important, 
kecause the profiles thus derived were later to be used as 
inputs tc low frequency (<1000HZ) acoustic models and the 
essential shape cf the profile was maintained. 
Additionally, ray theory is only valid if sound velocity 
does not change much in a wavelength, therefore, the inclu- 
Sion of the higher wavenumber fluctuations of temperature 
(and, therefore, sound speed) with depth could invalidate 
the use of this theory [Urick, 1983]. 

Four features of the vertical temperature profile were 
selected as being of primary importance in describing the 
thermal structure for low freguency passive acoustic appli- 
cations. These were: sea surface temperature (SST), mixed 
layer depth (MLD) (greviously defined), the below layer 
gradient (BLG) (defined as the temperature gradient in the 
25m immediately below the MLD) and the thermocline gradient 
(THG) (Defined as the temperature gradient between the MLD 
and the deep sound channel axis depth). The contour fields 
derived for each of these parameters, Figs. 4.7 to 4.10 , 


are descrited in detail in Section B of this chapter. 


Be. DESCRIPTION OF AREA THERMAL STRUCTURE 


The following is a brief description of the thernal 
structures measured in the study domain during the three 
CPTOMA5 cruises. The desciption in limited to the four 
parameters defined in section A above. 
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_ Figure 4.9 BLG Contours . 
a. Cruise A b. Cruise AII c. Cruise AIII. 
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Figure 4.10 IHG Contours. 
a. Cruise A b. Cruise AII c. Cruise AIII. 
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The significant weather for each cruise (i.e. wind speed 
and clcud cover) were obtained from R/V ACANIA's hourly log. 
The wind was from a generally southeastward throughout all 
three cruises, and, for the most part, wind speeds were 
between 15 and 20 knots. The major excursions from this norm 
were (a) at the beginning of cruise AI when wind speeds were 
greater than 30 knots on the first and second days (16 and 
17 Jun), (b) during the latter part of cruise AI and early 
in cruise AII when winds of 10 to 15 knots were recorded and 
(c)during the first two days of cruise AIII when wind sreeds 
again increased to 30 knots or more. Throughout the three 
cruises, partly clcudy to clear skies prevailed, with 
predominantly less than 4 oktas of cloud cover. Pressure 
remained nearly steady at approximately 1030 mb. 


1. Cruise AI (15 Jun 


SST*s ranged from 14.6°9°C in the northwest of the 
domain te 12.29°C associated with a cool feature pretruding 
into the southeast, Fig. 4.7 (a). The effect of tke cool 
intrusion was to intensify the horizontal temperature 
gradient in its vicinity and to, locally, distort the other- 
wise nearly meridional SST structure. The maximur temypera- 
ture gradient intensification occurred on the northern flank 
of the cool feature, where a gradient of 7.0°9C/100m was 
observed. The position of the cooler water agreed well with 
the position obtained from satellite imagery, Fig 3.6. 

The depth of the mixed layer, Fig. 4.8 (a) , varied 
from 20 to 55m ina complex manner throughout the area with 
the shallowest depths ina 50km-wide swath southeast to 
northwest across the centre of the area. No discernable 
correlation was observed between the SST and MLD patterns. 

Below layer gradients ranged from 3.59C/100m in the 
sector affected by the cold feature in the southeast to 


5-99C7100m in an area apparently associated with the 


OZ 


Shallower MLD's, Fig 4.9 (a). MThermocline gradients varied 
in a né€ar meridional manner with the weaker gradients in the 
west. The magnitude of the gradient varied from 1.5°C,/100m 
to 1.8°9°C,/100n, Fig 4.10 (a). 


2. Cruise AIT (29 Jun-4 Jul, 1983) 


A general warming trend waS apparent between the end 
of cruise AI and the beginning of cruise AII with SS1I's 
generally 2.0°9C warmer throughout the domain. There was also 
a distinct change in the SST pattern, Fig 4.7 (b). The cool 
feature initially cbhserved in the southeast sector had 
relaxed southwards and a second cool feature was observed to 
extend northward into the southwest of the domain. The 
warmest water remained in the far northwest sector, witha 
near meridional structure throughout the remainder of the 
area. Temperatures ranged from 14.4°9C to 16.29C. The warming 
trend was to be expected for the period involved. 
Climatclogy [Podeszoa, 1976] shows the SST increasing fron 
14.19C in June to 16.5°C in July. The mainly light to fresh 
winds and clear skies prevailing through the latter part of 
AI and early days of AII could be expected to produce the 
warming trend noted. Satellite imagery once more agreed 
well, gualititavely, with the position of the major surface 
features shown by the analysis, Fig 3.7. 

Associated with the warming trend and generally 
lighter winds, the MLD shoaled over the period retween AI 
and AII. Layer depths, some 20m Shallower than than those 
observed during AI, were measured with depths varying from 
near surface inthe northwest to 40m in the southeast but 
with a ccmplex structure throughout the area. 

The strengthening of the BLG, Fig 4.9 (b) compared 
to those observed during cruise AI indicated that the 
warming trend had been limited to the upper layers of the 
ocean and was, therefore, probably atmospherically driven. 
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Magnitudes ranged frcem 10.09C/100m in the west in asSsocia- 
tion with the lower SST's to 3.59C/100m in the centre of the 
domain, loosely associated with the deeper MLD's. An 
increase in THGts to ketween 1.79C/100m and 2.19C/100m with 
a generally similar pattern to that of cruise AI, emphasized 
the supposition that the warming had been limited to the 


upper layers of the ocean. 
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buise AIII (13) Jul- 139 only 1383) 


The SST analysis, Fig. 4.7 (c) showed the presence 
of cool water along the southern edge of the domain with the 
cool feature in the southwest during cruise AII, having 
relaxed southward and a further intrusion of cocl water 
extending northward into the area from the southeast. 
Satellite imagery agreed well with the position of these 
features, Fig. 3.8 .~. However, coverage was too intermittent 
to allow the feature inthe southeast to be positively 
related to the cool feature measured in a Similar position 
during cruise AI. Temperatures ranged from 12.4 to 16.0°9C 
with the warmest water remaining as an almost staticnary 
feature in the far northwest of the domain. A relatively 
intense horizontal temperature gradient of 10.0°C/100m was 
measured in a 30km-wide band around the cool feature in the 
southeast. 

MLD‘'s, Fig. 4.8 (c) were, on average, deeper than 
those measured during cruise AII, indicative of the higher 
winds prevailing during the early part of the cruise 
producing significant mixing in the surface layers. However, 
Shoaling of the MLD, in both the north and south of the 
area, tc 10m or less would tend to indicate that the wind 
mixing was restricted to the area being sampled by ACANIA 
during that period. The general horizontal structure of the 


MLD remained complex. 
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The magnitude of neither the BLG nor the THG changed 
Much from those measured during cruise AII. However, the 
distribution of the raximum and minimum BLG was’ signifi- 
cantly different from that of AII, although there was a weak 
increase in BLG in association with the cool feature in the 


southeast of the domain. 
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For each of the parameters discussed in the Charter IV, 
nanely, SST, MLD, BLG and THG, a difference field was 
produced by subtracting the magnitude of the parameter 
derived from the model analysis from that derived from the 
OPTOMAS data. The procedure was carried out for each model 
analysis at each of the 49 grid points in the domain. The 
difference fields thus calculated were contoured, Figs. 5.1 
Gna. 242 

Statistics were calculated from the contoured difference 
fields by examining the percentage of the total domain 
enclosed by discrete contour intervals in each difference 
field, and for each parameter. The *bin-size' for each field 
was established by examining the individual plots and 
subjectively selecting an appropriate difference interval 
such that each bin would represent between 10 and 20% of the 
dynamic range of the variable, whilst keeping the area 
covered by each bin large enough to make measurement frac- 
tical. The bin-sizes selected were: 

SST - 0.49°9C 

MLD - 10.0nm 

BLG - 1.0°9C/100m 

THG ~- 0.19C/100m 
The areas within each bin were measured using a Planimeter. 
Figs. 5.25 to 5.48 show the histograms of the distributicns 
obtained by the above method and Tables II to IV list the 
relevant statistics for each histogran. 

A negative value in the difference calculations implies 
that the model analysis was overestimating that particular 
parameter (i.e., Difference = Sea truth - Model value). 
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Because of a corruption in the data supplied by FNOC it 
was not possible to derive contour fields of the model anal- 
yses for the termination of cruise AIII. 

The larger grid spacing used in the TEOTS, hemispherical 
analyses (ca. 320kmn) meant that TEOTS comparisons were, 
necessarily, general in nature and should be viewed in terms 
of average values , since it i1s unlikely that such a large 
grid spacing will allow the analysis to retrieve the mesos- 
cale variability present in the area ({ horizontal vari- 
ability on a scale of the order of tens of kilometers). 

The following sections discuss the statistics obtained 
and qualitatively estimate the effectiveness of the model 
analysis in representing the sea-truth. A summary of the 
model performances is given in Section A.6 of this chapter. 
Throughout the Chapter, sea truth refers to OPTOMAS derived 
data. 


A. CCMPARISONS 


1. Cruise AI 


The analysis dates for cruise AI were 16 June, 1983 
(initial) and 23 June, 1983 (final). FNOC data files showed 
inputs into the the analyses for the study domain of three 
Bees LOr the initial, and twenty-three ET's (from R/V 
ACANIA) for the final analyses. The statistics for the 
difference fields derived for this cruise are given in Table 
II and the histograms from which these statistics were 
derived are shown in Figs. 5.25 to 5.32 


ae EOTS Analysis 


The initial SST analysis, Fig. 5.1 (b) showed a 
nearly meridional structure with the warmest water in the 
west of the domain anda weakly strengthening horizontal 
temperature gradient in the east. The difference field, Fig. 
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5.1 (c) showed the analysis to be overestimating the meas- 
ured SST for the majority of the area. The mean difference 
was -0.89C and the maximum difference was -2.09C associated 
with the intrusion cf cool water into the southeastern 
sector, Fig. 5.1 (a). The final analysis, Fig. 5.1 (d) 
showed a general movement of the analysis field towards the 
sea truth in both share and magnitude. The shape field was 
only marginaly improved and continued to neglect the pres- 
ence of the cool feature in the southeast. The mean differ- 
ence reduced to 0.1°C and the spread of the difference field 
about this mean was also reduced compared to that of the 
initial analysis. The maximum difference remained in the 
vicinity of the cool feature but with a reduced magnitude of 
-0.89. 

The MLD was overestimated in both initial and 
final analyses, Fig 5.2 (b) and (d). The slight shoaling of 
the layer between the two analyses (2m) had no significant 
effect on the difference fields, Fig. 5.2 (c) and (e). The 
modeled MLD was, on average, 18 to 19.5m greater than the 
sea truth and the maximum difference was between 30 and 40m. 
The near uniform structure of the modeled MLD in both anal- 
yses waS a gross misrepresentation of the complex structure 
evident in the sea truth contour field, Fig. 5.2 (a). 

The average magnitude of both the BLG, Fig. 5.3 
and THG, Fig 5.4 model analysis fields also moved closer to 
the sea truth between the initial and final model analyses. 
The BLG analyses, however, Showed a marked difference 
between modeled and sea truth Shape fields for both 


analyses. 
rk. TEOTS Analysis 


The initial TEOTS SSI analysis, Fig. 5.5, (3 
overestimated the sea truth in a Similar manner to that 


cbserved ii the equivalent EOTS analysis and by 
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approximately the Same average amount i.e., a mean 
difference of -0.9°9C. The mean difference was almost halved 
in the final analysis, Fig. 5.5 (c).~ However, the error in 
the modeled shape field remained constant as, therefore, did 
the difference distribution about the mean. The difference 
fields are shown in Fig. 5.5 (d) and (e). 

The modeled MLD was unresponsive to input data 
and showed a lack of features. The TEOTS analyses overesti- 
mated the MLD by an average of 12m. TEOTS MLD's were closer 
to the mean sea truth MLD than were the EOTS analyses but, 
rFecause of the lack cf Shape agreement between model and sea 
truth the distributicn of the difference field abcut the 
mean was Similar to that observed in the EOTS analyses. 

The BLG, Fig. 5.7 and THG, Fig. 5.8 appear to 
respond to the input data by moving closer to the sea truth 
in the mean. However, the lack of shape agreement between 
model and sea truth contour fields remained unchanged. The 
TEOTS analysis was initially further from the sea truth and 
responded more slowly to input data than the EOTS analysis. 


2. Cruise Ail 


The analysis dates for cruise AII were 30 June, 1983 
(initial) and 7 July, 1983 (final). The inputs to the anal- 
yses were 3 BT's and nil for the initial and final analyses, 
respectively. Statistics for the model/sea truth comparisons 
are given in Table III and the histograms from which the 
Statistics were derived are shown in Figs. 5.33 to 5.40 


a. EOTS Analysis 


The initial SST analysis, Fig. 5.9 (tr) showed 
Similar shape characteristics to those observed in the 23 
June EOTS analysis (cruise AT) and a general warming trend 
of 0.5°C was indicated. Over the same period the sea truth 
analysis, Fig. 5.9 (a) indicated an increase in SST of 2.0°C 
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and a marked change in the shape of the contour field. The 
initial difference field, Fig. 5.9 (c) showed the analysis 
to be an underestimate of the sea truth by an average of 
1.09C with a maximum difference of 1.6°9C and a general lack 
of agreement between model and sea truth shape fields was 
apparent. 

The final model analysis, Fig. 5.9 (d), showed a 
trend towards a meridional SST structure (away from the 
zonal structure apparent in the sea truth field) with no 
change in the average SST. The mean difference remained 
Similar to that calculated for the initial analysis. 
However, because of the increased disagreement between the 
model and sea truth shape fields, the spread of the differ- 
ence distribution about the mean increased. 

The FNOC climatological depth/temperature 
profile for the JuneyJuly period in the region of the study 
domain predicts a SST of 14.2°C which was in close agreement 
with the mean SST of 14.4°C observed in the model analysis. 

The initial MLD analysis, Fig. 5.10 (b) showed 
no Significant change, in Shape or magnitude from that 
observed in the 23 June analysis, Fig. 5.2 (d). However, 
changes in the sea truth field, Fig. 5.10 (a) e between 
cruises AI and AII meant that the difference field, Fig. 
5-10 (c), was markedly different from that of the 23 June 
analysis. Although the final MLD analysis, Fig 5.10 (d) did 
indicate a slight shoaling (5m to 10m), the MLD was overes- 
timated by the model analyses throughout. 

The shoaling of the MLD between the initial and 
final analyses noted above is compatible with the model 
analysis moving towards the FNOC July climatology for the 
region which shows a change in MLD from 55m in June to 37n 
in July. 

The model BLG's were overestimated and the shape 
of the ccntour field in poor agreement with that of the sea 
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EEuth. A smoothing of the modeled BLG field toa near 
constant value in the final analysis did, however, reduce 
the mean difference from the 1.29C/100m observed in the 
initial analysis to 0.39C/100m. The results of the compari- 
sons are shown in Fig. 5.11. 

Model and sea truth THG's were in good agreement 
for beth analyses, Fig. 5.12. 


bk. TEOTS Analysis 


The warming trend noted in the EOTS analyses was 
also ckserved in the TEOTS analyses, Fig. 5.13 The mean 
difference decreased from 0.4°9°C in the initial TEOTS anal- 
ysis tc near zero inthe final analysis. The horizcntal 
structure of the TECIS analyses remained meridional witha 
weakly strengthening horizontal temperature gradient in the 
east of the domain ccmpared to the complex near zonal struc- 
ture cbserved in the sea truth SST. 

Modeled MLD's, Fig. 5.14 remained unchanged by 
any Significant amount in both magnitude and structure fron 
those of the TEOTS analyses of cruise AI. The modeled MLD's 
were overestimated Ty between 5m and 35m, when compared to 
the sea truth. 

The BLG, Fig. 5.15 was overestimated by the 
model in both the initial and final analyses with the 
magnitude of the mean difference increasing by the final 
analysis. 

The THG, Fig. 5.16 was in reasonably good agree- 
ment with the sea truth throughout. 


3. Cruise AIII 


As previously stated, model analyses were not avail- 
able for the termination of this cruise. There were no BT 
inputs to either model for the 10 July, 1983 analysis. 
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The only significant change in either model was a 
general increase in the SST (0.5°C) indicated in the EOTS 
analysis. 

The significant changes in the difference fields, 
Figs. 5.17 to 5.24 were due to changes in the sea truth 
fields between cruises AII and AIII. 

The statistics of the difference fields and the 
relevant histograms are shown in Table IV and Figs. 5.41 to 
5.48 AS with cruise AII, the TEOTS SST and MLD analyses 
were, on average, in closer agreement with the sea truth 
than the FOTS analyses. 


4. 


To estimate the accuracy of the EOTS analysis in 
responding to input data the author's analysis program was 
revised to operate ona 4x4 grid with a 40km grid spacing, 
Similar to that used in the EOTS analysis progran. The 
twenty-three BT's input to the EOTS 23 June analysis were 
used aS input. The resulting SST contour field and that 
produced by the eguivalent EOTS analysis are shown in Fig. 
Sse 

From Fig. 5.49 it can be seen that the EOTS analysis 
did nct respond accurately to the input data. The accuracy 
of the author's analysis program had previously been veri- 
fied ty comparisons with manually contoured SST fields with 
favourable results. The reason for the lack of resolution in 
the EOTS analysis is not fully understood but is prokabkly 
due to excess weighting on the first guess field used in the 
FOTS FIB analysis. This first guess field provides the shape 
information to which any input data is fitted in a least 


squares sense. 
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Lack of synopticity in the OPTOMA5 data proved to be 
unimportant since the only inputs, of any Significance, to 
the FNOC models were the twenty-seven BT submissions frior 
to the 23 June analyses. The fact that the BT's were input 
to the models as one data set made them synoptic to the 
computer model and, therefore, a synoptic assumption is 
valid for the sea truth analysis with which the model anal- 


yses were compared. 
6. Summary 
ae OD 


The EOTS SST analysis responded to the magnitude 
Changes dictated by input data (when sufficient data was 
available) but horizcntal structure was not resolved ona 
scale commensurate with the resolution capability of the 
“Okm grid spacing used. In the absence of input data the 
EOTS SST moved towards climatology and became even less 
structured. 

There waS a measureable response to input data 
in the TEFOTS SST analysis. The magnitude of the response was 
less than that observed in the EOTS analysis and the lack of 
agreement in the share field between model and sea truth was 
even more pronounced, as might have been expected frecm the 
larger grid spacing used in the TEOTS analysis. ECL S 
perfomed better, in a general sense, than EOTS when data 
were not available and the move towards climatology noted in 
the ECTS analysis was not observed. 


be SiLD 


Neither model produced a realistic analysis of 
the sea truth MLD whether input data were available or not. 


Both analyses tended to overestimate the magnitude of the 


6 3 


MLD, with the 


estimate of the 


Cc. BLG 
The 
truth MLD made 


parameters were 


agreement noted 


TEOTS analysis giving a 
sea truth MLD than the EOTS analysis. 


consistently better 


and THG 


lack 
comparisons of 
defined 


between these two parameters 


and sea 
BLG and THG difficult (Both 
in terms of the Any 


of agreement between modeled 


layer depth). 
was therefore 


regarded as purely coincidental. 


B. INTER~ANNUAL VARIABILITY 

The mesoscale variability of the study domain has been 
dealt with in some detail in the previous sections. 
Underlying this mesoscale variability 1s a variability with 
a time scale of the order of years (the effect of El Nino). 
The region was known to be influenced by such an event 
during the period of the present study. To illustrate the 


effect of such an 
were compared for July, 


study domain. 


July climatology for the region. 


Was 


cruise of July/August, 


this cruise 


series of cruises) 


The mean profiles 


derived from 


was very close to 


anomaly, 
1983 


mean depth/temperature frofiles 
1982 in the 
were also compared to FNOC 


and July/August, 


The mean profile for 1982 
the OPTOMA2 
(The area of the CCS sampled by 
the OPTOMAS 


derived from XBT 


data 
1982 


XBI collected during 


that sampled in 
and that for 1983 was 


data collected during cruise AIII of the OPTOMA5S' series 
(July, 1983). Figs. 5.50 to 5.52 Show the results of the 
comparisons. 

FNOC climatology was in reasonable agreement with the 
1982 mean profile {assumed non-El Nino year), Fig 5.50 . the 
major discrepancy between climatological and measured 
profiles occurred in the upper part of the water column 


(surface to 50m) 


which could be expected since this is the 
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region subject to maximum atmospheric forcing and, there- 
fore, the short time scale variability impossible to predict 
with climatology. 

The presence of a thermal anomaly in the depth/ 
temperature profile was clearly evident when the _ 1982 
profile was compared to that from 1983, Fig. S220.) ihe 
assumed El Nino event was present asa layer of warmer 
waters between the near surface and approximately 200m with 
a Maximum anomaly of 2.59C at a depth of 60m. The inability 
of a static climatology to represent such an event was 
demonstrated by comparing the 1983 mean depth/temperature 
profile with that obtained from FNOC climatology, Fig. 5.52. 

To oktain a measure of the operational significance of 
an ancmalous event such as the one described above transnis- 
sion loss curves were computed from the FACT9H transmission 
loss model for theJuly 1983 climatological and mean depth/ 
temperature profiles . The only measurable effect, Figs. 
5.53 to 5.55 was a reduction i1n convergence zone range in 
the curves derived frem the climatological profiles (apfrox- 
imately 3%) and a decrease in convergence zone gain 


(approximately 5db). 


C. RESULTS OF RECENT FNOC STUDIES 


Recent modifications have been made to both EOTS and 
TEOTS analysis models (Frost, pers. comm.) and, in partic- 
ular, the model sensitivity to input changes in PLD (MLD) 
has been modified. The incorporation of the TEOTS analysis 
into the finer grid regional model has also been achieved 
for certain areas. The effect of these modifications is the 
subject of an ongoing investigation by FNOC, based on inputs 
from the Bay of Biscay and the Western Mediterranean Sea 


regicnal analyses. 
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The initial results from the study are promising with 
TEOTS apparently showing considerable skill in predicting 
Changes in layer derth even in the absence of input data 
(Frost, fers. comm.). No information is available from this 
study related to the skill of the analysis in resolving the 
horizcntal structure in the region. 
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Figure 5.2 Contoured MLD Fields for EOTS/Cruise AI 
(contour units am). 
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Figure 5.3 Contcured BLG Fields for EOTS/Cruise AI 
(contour units 9C/1000m). 
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Figure 5.4 Contoured THG Pields for EOTS/Cruise AI 
(contour units 9C/1000m). 
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Figure 5.6 Contcured MLD Fields for TEOTS/Cruise AI 
(contour units m). 
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Figure 5.7 Contcured BLG Fields for TEOTS/Cruise AI 
(contour units °C/1000m). 
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Figure 5.8 Contoured THG Fields for TEOTS/Cruise ATI 
(contour units °C/1000m). 
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Figure 5.9 Contoured SST Fields for EOTS/Cruise AII 
(contour units 9Cx10). 
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Figure 5.10 Contoured MLD Fields for EOTS/Cruise AII 
(contour units nm). 
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Figure 5.11 Contoured BLG Fields for EOTS/Cruise AII 
(contour units 9C/1000m). 
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igure 5.12 Contoured THG Fields for EOTS/Cruise AII 
— (contour units °C/1000m).- 
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migure 5.13 Contoured SST PFields for TEOTS/Cruise AII 
(contour units 9Cx10). 
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Figure 5.14 Contoured MLD Fields for TEOTS/Cruise AII 
(contour units am). 
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Pigure 5.15 Contoured BLG Fields for TEOTS/Cruise AII 
(contour units 9C/1000). 
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Figure 5.16 Contoured THG Fields for TEOTS/Cruise AII 
(contour units 9C/1000). 
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Figure 5.17 Contoured SST Fields for EOTS/Cruise AIII 
(contour units °Cx10). 
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Figure 5.18 Contoured SLD Fields for EOTS/Cruise AIII 
(contour units m). 
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Figure 5.20 Contoured THG Fields for EOTS/Cruise AIII 
(contour units 9C/1000). 
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Pigure 5.21 Contoured SST Pields for TEOTS/Cruise AIII 
(contour units °Cx10). 
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Pigure 5.22 Contoured MLD Fields for TEOTS/Cruise AIII 
(contour units om). 


88 





= 
|| 
" 





0 et a , SE REE RD NAR SS RE ~ a 





® 100 
a. Cruise AIII Analysis 
| 





Q 100 


| b. TEOTS Analysis 10 Jul c. Difference Field 10 Jul 


a. 


Pigure 5.23 Contoured BLG Fields for TEOTS/Cruise AIII 
(contour units 9C/1000). 
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Figure 5.25 SST Dafference Distribution EOTS/Cruise AI 
ae 16 Jun Analysis b. 23 Jun Analysis. 
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Pigure 5.26 MLD Difference Distribution EOTS/Cruise AI 
a. 16 Jun Analysis b. 23 Jun Analysis. 
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ae 16 Jun Analysis b. 23 Jun Analysis. 
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Figure 5.33 SST Difference Distribution EOTS/Cruise AII 
ae 30 Jun Analysis b. 7 Jul Analysis. 
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ae 30 Jun Analysis b. 7 Jul Analysis. 
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Figure 5.40 THG Difference Distribution TEOTS/Cruise AII 
as 30 Jun Analysis b. 7 Jul Analysis. 
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Figure 5.41 SST Difference Distribution EOTS/Cruise AII 
10 Jul Analysis. 
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Figure 5.42 MLD Difference Distribution EOTS/Cruise AIT 
10 Jul Analysis. 
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Figure 5.43 BLG Difference Distribution EOTS/Cruise AII 
10 Jul Analysis. 
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Figure 5.44 THG Difference Distribution EOTS/Cruise AII 
10 Jul Analysis. 
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Figure 5.46 MLD Difference Distribution TEOTS/Cruise AII 
10 Jul Analysis. 
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Figure 5.47 BLG Difference Distribution TEOTS/Cruise AII 
10 Jul Analysis. 
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Figure 5.48 THG Difference Distribution TEOTS/Cruise AIlI 
10 Jul Analysis. 
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b) Author's Contouring Progran. 
(same inputs as EOTS 23 Juln Analysis). 
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Figure 5.50 Temperature Anomaly vs. Depth 
Jul 1982-FNOC Climatology. 
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Jul 


983-Jul 1982. 


117 


TEMPERATURE ANOMALY (CCSS-CLIM) 


TEMPERATURE ANOM. 
0.5 1.0 1.5 


(DEG C) 
0.0 2.0 





ee ee gE A A SS 


Figure 5.52 Temperature Anomaly vs. Depth 
Jul 1983-PNOC Climatology. 
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Figure 5.53 Transmission Loss Curves~-S/S Disposition 
ae Climatology b. OPTOMA5S AIII (mean profile). 
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Figure 5.584 Transmission Loss Curves-S/D Disposition 
ae Climatology b. OPTOMAS AIII (mean profile). 
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Figure 5.55 Transmission Loss Curves-D/D Disposition 
a. Climatology b. OPTOMA5S AIII (mean profile). 
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Comparison Statistics Cruise AITI 
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TABLE IV 


Comparison Statistics Cruise AIII 


MODEL EOTS TEOTS 
PARAMETER 
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std. deviation 16.90 


= 


© © 
= NO 
© >) 


}—= 
WW 
>) 


—" 
Ud 
— 
© 


> 


Nn 
© 
© 


skewness -0.70 
mean 0205 

BLG std. deviation I.33 1 be 

Coma median 0.50 

skewness 0.20 

mean a. 20 0.04 

THG std. deviation 0.14 0.14 
Ceo) median 0.20 

skewness -0.20 Lo 0,30 





124 


Aw. DEFINITION OF ACOUSTIC TERAS 


e Convergence Zone (CZ) - Region of focussing or convergence 
cf sound. 

e Direct Range - The shortest distance from the receiver at 
which the transmission loss exceeds the Figure of Merit. 

e Figure of Merit (FOM) - The maximum allowable one way 
transmission loss (passive sonars) for which a 50% prob- 
ability of detection is possible for a particular sonar. 

e Leakage - Loss of energy from a surface duct to the main 
sound channel. 

e Low-Frequency Cut-Off - Frequency below which no trapped 
normal mode may exist for a given channel or duct. 

e Sound-Speed Profile (SSP) - Depth/sound-speed combinations 
defining the vertical variation in sound speed. 

e Transmission Loss (TL) - Acoustic parameter which quanti- 
tatively describes the weakening of sound intensity between 
a point Im from the source and a point at a distance in the 
ocean. 

e Vertical Beam-Width - Included angle between the maxinun 
upward and downward fropagating rays. 


Be. ASSUMPTIONS 


The cnly variable in the range calculations was to be 
the SSP. To achieve this goal certain representative, 
acoustic parameters were assumed constant throughout and 
were used as standard inputs to all transmission loss 
calculations. 


e The ocean bottom was assumed flat. 
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e A standard bottom-lcss/grazing angle curve applied to all 

frequencies. 

e A flat sea surface (sea state 1) input to the FACT nodels 

(assumed flat in the FE model) 

e The FOM'ts used in the range calculations were &80db and 

90db. These figures are representative of the range of FOM's 

encountered in operational systems. 

e Source/receiver disfositions used were; 
20m/20m-Sshallow/shallow (S/S) 
20m/200m-sShallow/deep (S/D) 
200m/200m-deep/deep (D/D) 

The S/S disposition was intended to give a _ measure of 
the surface duct acoustic propagation; S/D of cross-layer 
propagation and D/D of deep sound channel propagation. All 
three combinations can be utilized by presently operational, 


shipkcrne systems. 


C. DERIVATION OF SOURD-SPEED PROFILES 


Depth/temperature profiles were constructed for both the 
EOTS and OPTOMA5 analyses at the positions of maximum error, 
ketween the two analyses, for SST, MLD, BLG and THG. The 
position of the maximum error was obtained from the differ- 
ence field discussed in Chapter V and the depth/temperature 
profiles were generated by utilizing the contour fields 
derived as described in Chapter IV. 

SSP's were computed from the temperature profiles using 
Wilson's Equation. The depth/salinity information was 
obtained from the DE STEIGEUR, deep CTD data (Chapter ITI) 

Only the data relevant to cruises AI and AII were 
utilized. Long range prediction (CZ) comparisons were made 
ketween profiles based on the OPTOMA5S, EOTS and TEOTS anal- 
yses. Direct range comparisons were limited to profiles 
based on the OPTOMA5 and EOTS analyses only. 
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The SSP's derived for each model analysis (and all 
Subseguent acoustic calculations based on those SSP's) were 
designated according to the parameter, and analysis date eg. 
SST1 refers to a profiles derived at the position of maxinun 
difference between mcdeled and OPTOMAS SST for the initial 
model analysis of the particular OPTOMAS5S cruise being 
studied; MLD2 refers to the profile derived at the position 
of maximum difference in MLD‘'s for the final model analysis 


of the cruise, etc. 


De. TRANSMISSION LOSS MODEL COMPARISONS 


Differences in model physics made direct comparisons 
ketween FACTOH and the PE Models unwieldy, as did limita- 
tions in the model computer programs. One of the major 
differences was that the PE model has a vertical beamwidth 
limitation which causes a program termination if exceeded. 
The limitation is a function of the SSP and, for the 
profiles utilized in this study, the maximum vertical bean- 
width which could be maintained was 20°. No such limitation 
exists fer the FACT9H model which sums acoustic intensities 
from all forward directions, equivalent to a vertical tean- 
width of 180°. The effect of this beamwidth limitation is 
that the PE model will predict consistantly greater trans- 
mission loss for a given range than will the FACT9H model, 
assuming all other inputs to each model are identical. 
However, the effect of SSP changes relating to temperature 
profile differences between modeled and OPTOMAS ocean 
thermal structures were computed on a model-by-model kasis 
and the results compared, Table V. An example of the trans- 
mission loss curves computed from the two models for the 
same SSF, Fig. Gra alis Clearly shows the effect of the 
differing model physics. 
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Table V shows the Parameter Error (difference between 
OPTOMAS and EOTS vaiuves for the parameter, expressed as a 
percentage of the OPTCMA5 value); the Range(R) (direct range 
(km) computed from the transmission loss curve based on the 
OPTOMAS data); Range Error. ({E) (difference between direct 
range obtained utilizing OPTOMA5S data andthat utilizing 
EOTS analysis data, expressed as a percentage of R). A 
negative sign indicates that the model analysis had produced 
an overestimate of that particular parameter. The compari- 
sons shown in Table V are for a 600Hz frequency. 

The lack of agreement between the two transmission loss 
models was apparent with the PE model predicting signifi- 
cantly shorter ranges than the FACT9H model. The PE model 
was also far more sensitive to variations in the input SSP, 
particularly when both source and receiver were shallow, and 
showed a dependence on the SSP for the S/D and D/D disfosi- 
tions, which, as discussed in later sections, was not the 
case for the FACT9H model. Convergence zone ranges comruted 
by both FACT9H and the PE model were in good agreement. 
However, the gain in acoustic energy at the convergence 
zones waS approximately 5db lower in the PE model computa- 
tions than those of FACT9H. 

The remaining comparisons discussed in this chapter were 


based on the FACT9H transmission loss computations. 


E. SUMMARY OF ACOUSTIC COMPUTATIONS 


When comparing ranges predicted utilizing measured and 
modeled profiles, an error of less than 20% was assumed 
operationally insignificant. This assumption was Lelieved 
reasonable when the uncertainties associated with the FACT9H 
acoustic model accuracy (eg., the lack of attenuation due to 


absorpticn and volumetric scattering) are considered. 
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Comparisons of CZ ranges computed using SSP*'s 
derived from OPTOMA5, EOTS and TEOTS data, TableS VI to IX 
revealed a relative lack of sensitivity of the CZ range to 
thermal structure variations of the order of magnitude 
present between modeled and sea truth analyses. The average 
CZ range computed from sea truth data was 50km, varying 
between 44km and 62kn. Ranges computed utilizing SSP's 
derived from FNOC model analyses were within 10% cf those 
computed utilizing sea truth data. The maximum errors were 


observed when both source and receiver were shallow. 
2. Direct Range Comparisons 


Subjective analysis of the transmission loss curves 
computed for the long range comparisons discussed akove, 
combined with an understanding of ray-theory physics, led to 
the direct range ccmparisons being concentrated on _ the 
surface duct propagation paths (i.e., shallow source and 
receiver). A limited number of comparisons were made between 
direct ranges obtained from the transmission loss curves 
computed utilizing OFTOMA5 and those utilizing EOTS derived 
profiles for cross-layer and below layer propagation paths. 
The results from these comparisons, Tables X and XI ccnfirn 
the supposition that the FACT9H computed transmissicn losses 
for the S/D and D/D fropagation paths are relatively insen- 
Sltive to SSP variations of the order of magnitude present 
ketween profiles derived from OPTOMA5 data and those from 
the EOTS analyses. In the remainder of this discussion 
"range error’ will refer to the difference between ranges 
calculated from transmission loss curves computed from sea 
truth profiles and those computed from EOTS profiles. 

The maximum range error for the cross-layer propaga- 


tion rath was 18% ona range of 8km, with the majority of 
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the range errors within 10% of a sea truth range of 8km at 
the 50Hz frequency. However, a weak frequency dependence was 
observed and errors cf the order of 15% on a range of 8km 
were measured at 600HZ, associated with the lower FOM. The 
Maximum error also cccurred at this frequency, associated 
with the lower FOM, Table X . 

The range errors measured for the below layer propa- 
gation paths were generally negligible (<5%) and independent 
of freguency and FOM, Table XI . However, for profiles 
SSTiand SST2 errors of the order of 15% on a range of 18km 
for the lower FOM and 15% on a range of 77km for the higher 
FOM, were measured. The presence of a weak shallow sound 
channel in the sea truth SSP is believed to re the cause of 
these significantly higher range errors and would also 
explain the lenger direct ranges associated with the sea 
truth profile (i.e., 18km compared to the average of 8km for 
the remaining profiles) 

For the surface duct propagation paths (S/S), trans- 
mission loss curves were computed at frequencies of 50, 300 
600 and 1000Hz for profiles associated with cruise AI and at 
300, 600 and 1000Hz for those associated with cruise AII. 
The 50Hz frequency was removed from the comparisons in 
cruise AII, Table XIII, once the range errors were shown to 
be negligible at this freguency for the cruise AI comfari- 
sons, Takle XII . 

From the range comparisons, Table XII and XIII it 
can be seen that at the higher FOM and for frequencies Lelow 
1000HzZ, the transmission loss is insensitive to SSP varia- 
tions of the order of magnitude observed with a range error 
of less than 10% on a mean range of 55knm. 

For range calculations based on the lower FOM, 
Significant range errors were observed for frequencies of 
300Hz and above. The range error was also observed to be a 


function of freguency. Plots of range error versus parameter 
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error for the four parameters being studied revealed a 
strong dependence of the range error on the MLD error, Fig. 
6.2 . The curves shown in Fig. 6.2 were fitted to the points 
in an RMS manner and do not imply a functional relationship 
between MLD and range error but are intended to shew a 
general trend. The scatter of the points in Fig. SEs sloyeel 
cate the complex interaction between the various parameters. 
From the curves it can be seen that the range error 
increased with increaSing MLD error and frequency, for 
frequecieS above approximately 300Hz. The presence cf a 
non-zero range error for a zero MLD error was a_ further 
indication of the effect of the other parameters. Picts of 
range error versuS errors in the remaining three parameters 
were also generated, Figs. 6.3 to 6.5 These curves appeared 
to indicate a marked dependence of range error on each 
parameter. However, ylots of MLD error versus the errors in 
the remaining parameters, Figs. 6.6 to 6.8 revealed that the 
major range error dependence was indeed on MLD error; e.g. 
the trend in range error when compared to SST error is in 
the same general sense as the trend in MLD error when that 


too 1S ccmpared to SST error. 
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Figure 6.2 Plot of MLD Error vs. Range Error. 
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Figure 6.3 Plot of SST Error vs. Range Error. 
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Figure 6.7 Plot of BLG Error vs. ALD Error. 
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Figure 6.8 Plot of THG Error vs. MLD Error. 
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VII. SUMMARY AND CONCLUSIONS 

The grid spacing used in the FNOC ocean thermal analysis 
models (1.e., 40Okm for EOTS and 320km for TEOTS) made it 
extremely unlikely that they would be capable of recovering 
the mesoscale, horizcntal variability (of the order of 30 to 
4Okm) dcminating the study domain. 

Within the limitations of the grid spacing used, the 
FOTS SST analysis responded well to input data. However, the 
analysis methodology is believed to have smoothed the input 
field to the extent that significant detail was lost. 

The hemispherical TEOTS SST analysis responded to data 
input in magnitude but was totally insensitive to the hori- 
zontal shape information provided by the input data. This 
lack of response is believed to be a function of the coarse 
grid sfacing of the TEOTS analysis field. 

Both EOTS and TEOTS were insensitive to input mixed 
layer depth (MLD) information. The EOTS MLD analysis 
remained almost constant throughout the study period in both 
magnitude and shape. The TEOTS MLD analyses did change over 
the study period and were, on average, in closer agreement 
with the measured MLD magnitudes, but had an almost constant 
shape field. The lack of detail in both EOTS and TEOTS MLD 
analyses was in marked contrast to the complex and variable 
measured MLD shape fields. The improved performance of the 
TEOTS analysis over the EOTS analysis is believed to be the 
result of the influence of the TOPS element on the analysis. 

The response of the models to data input for depths 
below the MLD was inconclusive since the parameters tested 
(BLG and THG) were functions of the MLD and, thus, errors in 
the MLD analysis would propagate into these remaining 
fields. 
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The FNOC procedure of using a single climatological 
depth/temperature profile for areas such as the study domain 
was Shown to be inaccurate. The study domain was shown to 
vary on a Secular time scale (El Nino effect) with a 2 to 
39°C thermal anomaly, compared to FNOC climatology, measured 
at depths of 50 to 100m during the study period. The effect 
of a poor climatology would be reflected in the model anal- 
yses in the absence of input data. 

For acoustic frequencies below 1000Hz and based on 
results from the FACTI9H Transmission Loss model, the errors 
in the EOTS and TEOIS analyses were insignificant for (a) 
long range (convergence zone) predictions, (b) direct range 
predictions for sonar systems with a figure of merit of 90db 
or greater, and (c) direct range predictions for cross layer 
and below layer source/receiver dispositions. The errors in 
the model analyses were found to be significant for sonar 
systems with a figure of merit of 80db or less when both 
source and receiver were shallow. 

The errors in modeled mixed layer depth were found to 
have the greatest effect on accustic prediction errors, and 
the effect was found to be directly proportional to acoustic 
frequency. 

A comparison between the FACT9H and PE Transmission Loss 
models revealed significant differences between the outputs 
of the two models. 

The PE model appeared far more sensitive to errors in 
the temperature profile than did the FACT9H model. The 
direct range errors obtained from the PE model outputs, 
based on modeled and measured temperature profiles, showed 
no dependence on the figure of merit and less dependence on 
the source/receiver disposition than did those based on the 
FACT9H model outputs. 

The disparity between FACT9H and PE model sensitivities 
nade an absolute determination of the operational 
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Significance of errors in the EOTS and TEOTS analyses 
impossible. However, the maximum effect of depth/temperature 
profile differences on direct range predictions was 
Significant, whichever transmission loss model was’. used, 


when both source and receiver were shallow. 
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VIII. RECOMBENDATIONS 

The staff at FNOC are presently restructuring the ocean 
thermal analysis model algorithms to enable the models to 
respond more readily to input data, particularly the PLD 
segment cf the analysis. The physics of the TEOTS analysis 
are also being incorforated into the finer grid, regional 
analyses. The effectiveness of such changes should be tested 
in a Similar areal manner to that of the present study, and 
for many more areas of the world's ocean, to establish the 
operational value of such an analysis systen. The mininun 
data input reguired tc give a representative analysis under 
varying conditions should also be established, in order that 
assets may be deployed effectively and economically. 

For regions where anomalous events are known to occur, 
the feasability of utilizing a quasi-dynamic climatology, 
based on hindcasting, should be studied, i.e., climatology 
could te updated when a persistent anomalous event is known 
to be present. This would give a more representative anal- 
ysis, at depths below those influenced by the TOPS elements 
of the TEOTS analysis, in the absence of input data. 

Satellite IR imagery could be a powerful asset in data- 
Sparse areas of the cceans, providing a necessary input to 
ocean thermal analysis models. However, the results of the 
present study have shown that the relationship between SST 
and sub-surface features is not always an obvious one. 
Studies should be undertaken to relate surface features with 
sub-surface structure, either through direct correlations or 
modeling. Such relationships would also assist the planning 
of sampling strategies aimed at achieving the most econon- 
ical and efficient use of assets whilst providing sufficient 
input to the ocean thermal analysis model to achieve the 


desired accuracy. 
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The effect on the higher fregquency;i.e¢., greater than 
1Khz, active sonar performance of errors in the model anal- 
ysis fields should be examined. The loss of vertical struc- 
ture could prove important. In particular, areas where 
Significant vertical features are present; e.g-, Shallow 
sound channels, should be examined and the feasibility of an 
additional floating level in the analysis, to represent the 
depth of such features, should be addressed. 

The effect of mesoscale variations, on scales such as 
those observed in the present study, on the other parameters 
in the sonar equation should be examined; @e.g.-, ambient 
noise changes due to changes in the marine ecosystem in the 
region of the weak ocean fronts associated with upwelling 
regions. 

Operationally important areas are, in general, the more 
densely sampled regicns of the ocean. However, the present 
fixed regional boundaries utilized in the FNOC analysis 
models are only broadly representative of these operational 
areas and generally extend beyond the limits of prime 
interest. AS a consequence, the shape field used in the 
analysis is influenced by possibly less well sampled areas 
outside the region of interest with possible erroneous 
smoothing effects on the horizontal and vertical thermal 
structure that the analysis would otherwise produce. To 
remove such a possibility, the utilization of a flexible 
regional boundary, fitted to the area of interest , and with 
a grid spacing reprasentative of the horizontal variability 
of the region, should be studied. Such an approach as the 
one outlined above would also facilitate the incorporation 
of a more detailed, possibly dynamic, climatology into the 


analysis. 
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